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Abstract 



We consider the proton decay in supersymmetric models with a gravitino 
or axino lighter than the proton. This consideration leads to a stringent 
limit on the R parity and B violating Yukawa coupling of the superpotential 
operator UfD^D^ as \'{ l2 < 10 _15 (m 3 / 2 /eV) for a light gravitino, and \'{ 12 < 
lCr 15 (F a /10 10 GeV) for a light Dine-Fischler-Srednicki- Zhitnitskii axino. For 
hadronic axino, the constraint is weakened by the factor of 10 3 . 
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Typeset using REVIgX 



Proton stability strongly constrains the baryon (B) and lepton (L) number violating cou- 
plings. Since all known fermions lighter than the proton carry a nonzero lepton number, the 
couplings (or the combinations of couplings) relevant for the proton decay should conserve 
B — L. However if there is a lighter fermion which does not carry any lepton number, proton 
decay may be induced by a B violating but L conserving interaction alone |TJ]. There are in 
fact very interesting class of models which predict such a light fermion. In supersymmetric 
models in which supersymmetry (SUSY) breaking is mediated by gauge interactions, the 
squark and/or gaugino masses, i.e. the soft masses in the supersymmetric standard model 
(SSM) sector, are given by m so f t ~ {-) n Ag where n is a model-dependent positive integer 
and As corresponds to the scale of spontaneous SUSY breaking 0. In such models, in order 
for m soft to be of order the weak scale, A s is assumed to be 10 ~ 1000 TeV, leading to 
the gravitino mass m 3 / 2 = A s /Mp < 1 keV far below the proton mass. If a global U(1)pq 
symmetry is introduced in a gauge-mediated SUSY breaking model to solve the strong CP 
problem by the axion mechanism ||, SUSY breaking in the axion sector is mediated also 
by some gauge interactions. The axino mass in such models is given by m« ~ (a/7r) m A s /F a 
where m is again a model-dependent (but typically bigger than n) positive integer and F a 
denotes the scale of spontaneous U(1)pq breaking [f|. Obviously then the axino is lighter 
than the proton for a phenomenologically allowed F a > 10 10 GeV. In other type of models 
in which SUSY breaking is transmitted by supergravity interactions, the gravitino mass is 
fixed to be of order the weak scale, however there is still a room for an axino lighter than 
the proton ||. As was pointed out in Ref. ||, some supergravity- mediated models lead to 
m- a ~ fn 3 / 2 ( m 3/2/Mp) 1 ^ 2 — 1 keV for which the axino would be a good warm dark matter 
candidate 0. In this paper, we wish to examine the proton decay involving a light gravitino 
or axino to derive a constraint on the superpotential interaction X"j k UfD^D^ which violates 
R parity and B, while conserving L. 

Let us first consider the proton decay involving a light gravitino, more precisely the 
helicity ±1/2 Goldstino component. Our starting point is the effective lagrangian below the 
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scale As = Jfri^^Mp but above the weak scale soft mass m so f t : 

£ = £ssm + £g, (1) 

where £ssm denotes the lagrangian density of the SSM fields and the Goldstino lagrangian 
C G is given by M 



Cg = l&fdpG + ^ (\ a Ya^d p GF; u + 2^,(1 - TshV^GA,^ + h.c) (2) 



i 

where G denotes the four- component Majorana Goldstino field. Here £ssm includes the 
terms associated with the B violating superpotential interaction, 

V^ssm 3 X" jk U t c Dp c k , (3) 

and (<f>i,ipi) and (\ a ,F® u ) stand for the left-handed chiral matter and gauge multiplets in 
the SSM sector. Note that the above form of Goldstino lagrangian is enough for the study 
of the process involving a single on-shell Goldstino obeying i'j^d^G = m^iiG. 

Integrating out all fields heavier than the scale of the QCD chiral symmetry breaking, i.e. 
A x ~ 1 GeV, we are left with an effective lagrangian of the light quarks, q a (a = (u, d, s)), 
and gluons together with the light Goldstino (of course also the light leptons and the photon 
which are not relevant for our discussion). The operators responsible for the proton decay 
in this effective lagrangian at A x are induced by the exchange of the SU{2)l singlet squarks 
as 

O cS = ^T V % (5,(1 - riGWMl - l^G). (4) 
Here ml denotes the squark masses which are assumed to be (approximately) universal, 

Vdsu = Uuds = Uusd = 1) (5) 

and all other components of y a p^ do vanish. Note that the above operator has B — S — — 1, 
and thus the relevant proton decay mode is p — > G + K + . For a generic non-universal squark 
mass matrix, S = operator can be induced also to give rise to p — > G + tt + , however it is 



suppressed by a small squark mixing. To arrive at the above interaction operator, we have 
used the equation of motion of the on-shell Goldstino field and ignored the piece suppressed 
by the small 7713/2. Also ignored are the renormalization effects between the weak scale and 

Ax- 

The hadronic matrix elements of the above B = S = — 1 operator would be described by 
an effective chiral lagrangian including the Goldstino field. Let us consider a chiral operator 
O x which would induce p — > G + K + as a low energy realization of the light quark operator 
O e g below A x . Obviously it can be written as O x = Z^{\ — 'y^d^G where Z^ is a fermionic 
B — S — — 1 operator including P and K + . If Z^ does not include any spacetime derivative, 
O x is suppressed by the small factor m 3 / 2 /m p (for on-shell Goldstino) where m p denotes the 
proton mass. For containing a single spacetime derivative, we have 

O x = R 2 l Kl2 - (P(l - i,WG)d»K\ (6) 

where again the equations of motion are used together with "C m p . To estimate the 
size of the hadronic coefficient £, we use the naive dimensional analysis (NDA) rule of Ref. 
0, yielding 

\i\ * ±Kfl (7) 

where f n = 93 MeV is the pion decay constant. In fact, the NDA rule gives A x = 4irf n and 
then the typical energy in the proton decay, i.e. m p , is comparable to A x . This means that, 
within the NDA rule, chiral operators with more spacetime derivatives are equally important 
as the operator of Eq. (6). However for an order of magnitude estimate of the hadronic 
matrix element, the consideration of Z^ with a single derivative would be enough. Then 
applying the experimental limit on p —>■ K + + v for p — > K + + G induced by the interaction 
of Eq. (6), we find the following constraint on the R parity and B violating coupling: 



which is one of the main results of this paper. 



Let us now consider the proton decay involving a light axino. Similarly to the case of a 
light gravitino, we start from the effective lagrangian at scales below the scale F a of U(1)pq 
breaking but above m so ft'- 

C = Cssm + £a, (9) 

where the axino lagrangian Ca can be read off from 

J d 2 6d 2 6^(A + A^\<S>j + y (pe^pAW a W a + h.c}, (10) 

where A = (s + ia) + \p2Qd + F^O 2 is the axion superfield containing the axion a, the saxion 
s and the axino a, while <3>/ and W a are the chiral superfields for the SSM matter and 
gauge multiplets (</>/, ijji) and (\ a ,F® u ), respectively. Here c 1 and c a are dimensionless real 
coefficients. For F a defined as the scale of spontaneous U(1)pq breaking, the coefficients 
c a of the axion coupling to the gauge multiplets are of order unity in general. However 
as we will discuss later, the size of the coefficients c f of the axion coupling to the matter 
multiplets is somewhat model-dependent. Note that the above lagrangian corresponds to 
the super symmetric generalization of the conventional axion effective lagrangian 0: 

C a = ^d^n^j + ^^-aF^F^. (11) 

Obviously it is manifestly invariant under the nonlinear U(1)pq transformation, A — > A + ic 
(c = real constant), up to the PQ anomaly. At any rate, the relevant axino lagrangian is 
given by 

L A = l -ra^d^~a - -^-(^^(1 + 75 )S0* + h.c) 

l It a 

+ i5fe™ 1 -^ + Lc) ' (12) 

where a denotes the four-component Majorana axino field. Again the exchange of the SU(2) l 
singlet squarks leads to the following B = S = — 1 interaction in the effective lagrangian at 



O eS = lKl2 l a ^ {q a (l - 7 5 )^g 7 7"(l + 75)5, (13) 



where c 7 (7 = u, d, s) denotes the axino coupling to the supermultiplet containing the SU (2) l 
singlet right-handed light quark q lR in Eq. (12) and the squark degeneracy is assumed also. 

Similarly to the gravitino case, in order to estimate the proton decay rate from the above 
effective interaction, we consider a chiral operator of the form O x = X M 7 M (1 + 75)0 where 
X M is a B = S = — 1 fermionic current made of P and K + which are on mass-shell. For 

oc K + P r y^ 1 the chiral operator O x with the smallest number of spacetime derivatives is 
given by 

^^(P(l +75 )a)^, (14) 
where the hadronic coefficients £ 7 are again determined by the NDA rule as 

lerl * 16tt 2 ^. (15) 
This then leads to the experimental bound on the R parity and B violating coupling as 

^^o-.m'f^m, (16) 



.300 GeVj \ c 7 |£ 7 | J V10 10 GeV, 
which is another result of this paper. 

The above constraint from the proton decay involving a light axino depends upon the 
dimensionless coefficients c 7 describing the axino coupling to the supermultiplets of the 
SU(2)l singlet quarks [see Eq. (12)], as well as the axion scale F a . In fact, the size of c 7 
has a certain model-dependence. If the quark superfields carry a nonzero U(1)pq charge, 
which would be the case for the supersymmetric extension of the Dine-Fischler-Srednicki- 
Zhitnitskii (DFSZ) axion model [[Uj], the coefficients c 7 would be of order unity in general. 



However in hadronic axion models [11| in which all SSM fields have a vanishing U(1)pq 
charge, the coefficients c 7 are zero at tree level. However the axino-quark couplings are 
radiatively generated through the axino coupling to the gluon multiplet, yielding c 7 ~ 
(a c /n) 2 ln(F a /m so { t ) ~ 10~ 3 ~ 10~ 2 0. Thus the constraint for hadronic axion models 
becomes weaker than that for DFSZ models by the factor of 10 2 ~ 10 3 . 

To conclude, we have considered the proton decay involving a gravitino or axino lighter 
than the proton. Generic models in which supersymmetry breaking is mediated by gauge 

6 



interactions contain such a light gravitino. Then the R parity and B violating coupling 
A" 12 is strongly constrained by the proton stability [see Eq. (8)] to be less than about 
10~ 15 (m 3 / 2 /eV). About the possibility of a light axino, gauge-mediated supersymmetry 
breaking models endowed with a global U(1)pq symmetry generically predict an axino lighter 
than the proton. Also some supergravity-mediated models can give rise to a light axino, while 
the gravitino mass in such models is fixed to be the weak scale. We find that A" 12 in models 
with a light axino is constrained [see Eq. (16)] to be less than about 10 _15 (F a /10 10 GeV) 
and 10~ 12 (F a /10 10 GeV) for Dine-Fischler-Srednicki-Zhitnitskii axion models and hadronic 
axion models respectively. 
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